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Abstract 

' Gravitational field of stars of the Galaxy causes additional mo- 



tions of images of extragalactic sources. For typical realization of 
stars this results in a small rotation of extragalactic reference frame 
s * \ in the direction of Galactic rotation. We estimate angular velocity of 

corresponding image motion and show that it is different in case of 
<^ ■ discrete and continuous matter distribution in the Galaxy. 

q ! 1 Introduction 

Perspectives of microarcsecond astrometry raise questions concerning mo- 
tion of images of extragalactic radiation sources due to gravitational field of 
moving stars pQ , [2j , [Sj , [I] , jH] • This gravitational image motion (GIM) may be 
comparable to proper motion of quasars. There is some hope that correspond- 
ing accuracy to observe such effects will be available for future space-based 
radio interferometers pQ. However, it will be difficult to separate GIM from 
real proper motions, and this puts some limit for accuracy of fundamental 
reference frame based on extragalactic sources Indeed, we cannot take 

into account positions of all stars in the Galaxy; we only may work with prob- 
ability distribution of Galactic stars leading to distribution of image motions 
of certain extragalactic object 0, [Hi- 
lt has been pointed out in 0,0 that stochastic GIM of an extragalactic 
source induced by stellar motions is accompanied with a systematic compo- 
nent which depends upon bulk velocity of microlensing stars, including stars 
that are far away from the line of sight. This requires some explanation. If 
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we perform a statistical averaging using all possible positions of microlenses 
(stars), we shall obtain a value different from the most probable value of GIM 
that might be observed at the present epoch. There is also a difference be- 
tween the statistical average and arithmetic mean over all reference sources 
used in extragalactic reference frame. This difference appears to be essential 
(an example will be given in this paper below). The reason is that in reality 
the number of these sources is not sufficiently large. In this paper we treat 
this question in more detail. We show that one can estimate probable GIM 
effect, without recourse to probability distributions, after some restriction of 
sample of observed sources. The result differs in case of continuous (e.g., 
dark matter) and discrete (e.g., stars) mass distributions. 

2 Motion of source image microlensed by a 
point mass 

In case of microlensing of quasars by Galaxy stars we may neglect cosmo- 
logical curvature. Evidently this does not contradict to assumption that the 
radiation source is at the infinity. 

Let unperturbed light ray moves from the infinitely distant radiation 
source in negative direction of z-axis of Cartesian coordinates {x, y,z}, the 
observer being at the origin. Let position of microlensing point mass M be 
(r, z), where r = (x,y) is a two-dimensional vector in the transverse plane; 
i.e. r = |r| is the impact distance of the unperturbed ray with respect to the 
mass. We consider the case of weak microlensing r >> (mz) 1//2 , m = GM/c 2 ; 
therefore the light trajectory in the gravitational field of the point mass can 
be obtained in the post-Newtonian approximation. Note that in our problem 
we may have r of the same order as z. In this case 
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where = SE^) is a two dimensional vector describing the source image 
angular shift per unit value of m [H| , [I| , [I| • Formula (JTJ) may be obtained 
from a more general relation [7] in the case of an infinite source. 

Further V p = (v, w) stands for velocity of the point microlens, w is 
the velocity component parallel to the line of sight and v represents the 
transversal components. In virtue of (JTJ this leads to the source GIM that 
equals to mU (in radians per unit of time), where [3] 
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Here U is a function of the microlens position (r, z) and its velocity v. 

3 Motion of source image microlensed by Galaxy 
stars 

Because of smallness of the effect, the action of all Galaxy stars will be taken 
into account linearly. This will be performed by integrating (J2J). However 
this requires justification. In fact we are dealing with extended sources; this 
is not taken into account by Eq.(j2J). As we shall see below, this may be 
important in our problem. 

We are interested in some estimate of probable GIM value that may 
be obtained after taking arithmetic mean for as many reference sources as 
possible. As we pointed out in the Introduction, this is not a statistical 
average, because in reality the number of observed sources is limited. For 
example, the total number of extragalactic sources in ICRS is about 600. 
Moreover, this is much lesser, if we confine ourselves, e.g., to the Galaxy 
plane. 

We deal with an ensemble of possible positions of stars in the Galaxy. 
Consider a remote extended source, which has impact distance with respect 
to z-th star in this realization. Let p = minjrj}, where this minimum is 
taken over all Galactic stars. Now we consider two type of events: (A) when 
p » L, where L = D ■ as, as is angular size of the source, D is a typical 
distance to microlens (of the order of 50 kpc; (B) when p ~ L or p < L (at 
least one of stars is projected onto the source). We suppose that for typical 
value of L we have L < Re = (4ml)) 1 / 2 , in this case we may strengthen 
our arguments by considering even larger domain B', defined by condition 
p < L for corresponding realizations of stars. The rest of events we refer to 
the domain A'\ corresponding probability is Pa> ~ 1- The events from B' 
can be separated in observations: these are strong microlensing events that 
are characterized by considerable brightness amplification and relatively fast 
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image motions. These events might be taken into account, including the 
case of an extended source 0,0- Probability of these events Pb> ~10~ 6 is 
so small that they are not practically essential. In case of limited number of 
reference sources it is reasonable to exclude such rare microlens realizations 
that would hardly occur during this century. 

Nevertheless, though P B i « Pa>, the input of B' into the statistical 
average GIM < U > is essential because the events B' introduce large image 
velocities. To present an example, we derive below a nonzero value of GIM 
for a fictitious observer at the center of the Galaxy. Because of stationary 
rotation of the Galaxy, one must have < U >= (no matter, for an infinite 
time average or for a statistical average). This shows that in this case the 
sets A' and B' compensate each other. However, for a limited time interval 
we typically have a nonzero GIM value, because we do not meet realizations 
from B'. 

Therefore, the question is how to obtain a consistent GIM estimate. The 
most correct way is to obtain probability for all necessary velocity intervals. 
However, to have an order-of-magnitude estimate, we can avoid calculations 
of probability distributions, if we confine ourselves to obtain average value 
of GIM in the domain A'. Further we have in mind just this average when 
speaking about average GIM effect. It is important to note that this value 
practically does not depend upon the exact size of B'. This is provided by 
convergence of integrals that will be considered below. 

Thus we return to Eq.(j2J), which describes contribution of a single star. 
We assume that a star at the point (r, z) has velocity V p (r, z). Then we 
consider sum over all stars in the Galaxy, which must further be averaged 
with the Galactic mass density p(r, z). This enables us to pass on integration 
yielding average GIM in the domain A' for the weak microlensing events. 

< U tot > A >= ^ J dz J d 2 r p(r, z) U(r, z, V p (r, z)) (3) 

where we suppose that p vanishes outside a bounded domain. 

Taking into account the explicit form (j2J) it is easy to see that considerable 
contribution may be due to stars at large impact distances from the line of 
sight. Also one can show that the singularities in (jSJ) for small r are integrable. 
This allows us to avoid the question about exact value of lower limit of r in 
the definition of the domain A'. 
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4 Continuous versus discrete mass distribu- 
tion 



It is interesting to compare GIM (jSJ) with a corresponding expression in case 
of a continuous mass distribution p = p(r,z,t), satisfying the continuity 
equation 



dp 
dt 

In this case Eq. (0) should be integrated with the density p: 



-£ + div{pV p ) = 0. (4) 



G r r 
®tot = -2] dz j d 2 rp(r,z,t)^(r,z). 



Note that we have the same formula for ^tot in case of continuous and 
discrete matter distribution. However, this is not true for GIM. 

The shift of remote source image in case of continuous matter distribution 
has been treated in These shifts appear to be of the order of 10~ 5 4- 10~ 6 
radians. However this value cannot be observed from the Solar system. In 
order to deal with observable (in principle) values one must derive changes 
of the source positions with time: 

This will be compared to < XJtot >A' ■ After standard elimination of 
singular points at r = in U^ o4 and taking into account Eq. (JU) we represent 
this integral as sum of two terms 

UL = u<nu0 (6) 

being two-dimensional vectors, the components of the first term 

(i=l,2) 

Uf ] = ^JdzJd 2 r p(V p V) % 
are the same as the components of Eq. 

U<°> =< V tot > A , . (7) 

The second term is 
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U] = -— j dz J d 2 rdiv(V p p^i). 
For this latter term we have in cylindrical coordinates {r, ip, z} 
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= -§ J dzjd<pjdr{£ {rpV r Vi)} = % J dzj dip lira (rpV r %) , 
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(8) 

where V r = (V p • e r ), = (V p • e^), e r = {cos(y?), sin(<^), 0}, e v = 
{— sin(</?), cos (</?), 0}, is the longitudinal component of V p along z axis, 
the latter is directed along the line of sight to the source. For r^Owe have 
r\I> — ► — 4e r (2 > 0), and r^S? — > (z < 0). Then simple calculations yield 

. „ 00 



where represents transverse components of V p . Thus is the differ- 
ence of GIM effects for continuous and discrete distributions defined by mass 
density on the line of sight. Note that does not go to zero in case of mi- 
crolens masses fragmentation, if we make these masses smaller with leaving 
p unchanged. 



5 Average image motion in case of our Galaxy 

We use a four-component model of mass distribution in the Galaxy (buldge+ 
disk + halo + dark corona) according to [Bj. We confine ourselves to this 
model in view of its simplicity, though more recent and more adequate models 
of Galactic density are available at present. However, for order-of-magnitude 
estimates the model of |S| is enough. In this model the spherical components 
(bulge, halo and corona) have isothermal sphere mass density 

pM- 3MsRl 



4tt (R 2 s + R 2 ) 
R being distance from the Galactic center. 
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The parameters of spherical components are taken as follows: M B = 1.5 • 
10 11 (in Solar masses), R B = 5kpc; M H = 5-10 10 , R H = 25kpc; M c = 8-10 11 , 
Rc = 50kpc. We changed these parameters as compared to [S| to have better 
correspondence to the observed rotation curves. 

In case of the Galactic disk (in the cylindrical coordinates r, </?, z with the 
origin at the Galactic center) 
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where H = OMpc, R D = 15kpc and M D = 8 • 10 10 . 

We use this model along with corresponding rotation curves, which are 
recovered by means of potentials of the above mass distributions [Hj. The 
dependence of the rotation velocity around the Galactic center in the Galaxy 
plane is 



V 2 = G 



M B r 2 
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M c r 2 
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M H r 2 



(R 2 B + r 2 )^ (R 2 c + r 2 )^ (R 2 H + r 2 )i 2H 



n(r) = 2r ^ _ r ° + H r o~ H 

sjr 2 + r 2 D y/r 2 + (r D + Hf ^/r 2 + (r D - H) 2 

For an observer at the Galactic center we have U£ oi = (see Eq.(jHJ); this 
is evident in case of a stationary mass density and the observer at rest. We 
assume that all the components except continuous corona consist of stars 
(though if the corona also consists of stars, this does not change the re- 
sult significantly). Then in view of Eqs.(|B|l , (JJJ) we have the average GIM 
< \Jtot >A ,= — U*- 1 -*. This yields for source in the Galactic plane 3 • 10~ 8 arc- 
seconds per year. This also can be checked directly by means of the formula 
Q. When inclination increases, the effect decays strongly due to decreasing 
of star number on the line of sight. 

If the observer is located in the Solar system, we must take into account its 
own velocity and the velocity of stellar rotation around the Galactic center. 
Note that it is more convenient to calculate first U^ ot and , then the value 
of average GIM < U tot > A i is defined from Eq.tJHJ). Dependence of GIM upon 
the Galactic azimuth of sources in the Galactic plane for the observer in the 
Solar system is shown on FigHJ 
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Figure 1: Observer in the Solar system, sources near the Galactic plane: 
average GIM against Galactic azimuthal angle. 
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6 Discussion 



We have shown that GIM differs for a continuous and discrete mass density 
distribution. In the latter case GIM performs random walks [2],|II],|I] as 
distinct from a regular motion in case of continuous matter. In this paper we 
draw attention to the fact that the difference between these two cases remains 
also after averaging of the stochastic GIM (for the same mass density). For 
example, the observer at the Galactic center would see an additional fictitious 
motion of quasars around this center. However, if most of the Galactic matter 
would be continuous, the GIM effect will be absent. 

For an observer in the Solar system we have a nonzero average GIM 
value leading to apparent rotation of extragalactic reference frame in the 
direction of the Galaxy rotation. We stress that the averaging procedure 
involves only those realizations of Galactic stars, which have concern with 
weak microlensing. For sources in the Galactic plane maximal effect amounts 
about 1.5 • 10~ 7 arcseconds per year. This is extremely small. However, in 
principle, the effect can be observed, because it falls off with inclination of 
the line of sight to Galactic plane. Formal algorithm for such measurement 
must involve observations of sources at different Galactic latitudes. In this 
procedure the strongly microlensed images must be eliminated. 

As concerned practical measurement of the average GIM effect, we note 
that it is far beyond modern possibilities. Moreover, observation of the aver- 
age effect requires too many "good" extragalactic sources with approximately 
the same Galactic latitude. This is necessary in order to separate GIM from 
random image motions that may be much more essential. 

This work is supported by Science and Technology Center in Ukraine 
(Project NN43). 
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